Nonparaxial spatial optical solitons with nonlocality 

nonlinearity 



Jilong Chen, Zhiwei Shi 

1 Faculty of Information Engineering, Guangdong University of Technology, Guangzhou 

510006, P.R.China 

Huagang Li 

2 Department of Physics, Guangdong University of Education, Guangzhou 510303, 

P. R. China 



Abstract 

We investigate one- dimensional nonparaxial spatial optical solitons with non- 
local nonlinearity. We show an exact analytical solution to the nonlocal 
nonlinear nonparaxial propagation equation in the cases of high and weak 
nonlocality. We also numerically find that the degree of nonlocality can af- 
fect the width of nonlocal soliton beams, but have no effect on their stability. 
Contrarily, nonparaxiality can affect their stability, but have no effect on 
their width. 
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1. Introduction 

A nonlocal nonlinear response played an important role on the optical 
spatial solitons over the years. Nonlocality exists in different physical media, 
as nematic liquid crystals [lj, photorefractive media 0], thermal jsl, 0], and 
so on. Various degree of nonlocality given by the width of the nonlocal re- 
sponse function R(x) and the intensity profile of the beam can be divided into 
four types, like local, weakly nonlocal, general nonlocal and highly nonlocal 
response [5[]. 
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Nonparaxiality can refer to two different contexts of distinct nonparaxial 
character: high intensity and large angles of propagation p. The first type of 
nonparaxiality results from the evolution of ultra-narrow beams in nonlinear 
media. In 1993, Akhmediev et al cast doubt on the suitability and limitations 
of the normalized nonlinear Schrodinger (NLS) equation for describing the 
evolution of such beams jz|. The scalar theory of the self-focusing of an 
optical beam is not valid for a very narrow beam, and a vector nonparaxial 
theory is developed from the vector Maxwell equations [§, Based on 



these equations, Crosignani et al have reported and analysed bright [lOj and 
dark 11] nonparaxial solitons. In contrast to the first type, the second type 
of nonparaxiality arises from the rapid evolution of the field envelope of 
a broad (when compared to the wavelength) beam propagating at a large 
angle to the longitudinal axis [6j. The scalar nonlinear Helmholtz (NLH) 
equation well describe this nonparaxiality and overcome the limitations of 
Exact analytical soliton solutions have been found in a focusing 
Nonparaxial theory based on NLH equation has also 

boundary [16 



the NLS 12 



Kerr type medium 13 



been applied to find dark Kerr [14 



15 



and 



two- comp onent 

bistable |l7| Helmholtz soliton solutions. At nonlinear interfaces, Soliton 
refraction effects have a strong inherent angular character and constitute an 
excellent testbed for nonparaxial Helmholtz theory 

In this paper, we investigate the nonparaxial solitons in the nonlocal non- 
linear medium. We find an exact analytical solution to the nonlocal nonlinear 
nonparaxial propagation equation in the cases of high and weak nonlocality. 
In addition, we numerically study the influence of the degree of nonlocality 
and nonparaxiality on the width and stability of solitons. We devote Sec- 
tion [2] to a detail description of the theory model of the nonlocal Helmholtz 
soltions. Section [3] studies the exact analytical nonlocal Helmholtz solitons 
solutions, which correspond to the high and weak nonlocality. Section H] is 
devoted to studying the numerical solitons soluitions. The aim is to illustrate 
the influence of the degree of nonlocality, nonparaxiality on the width and 
stability of solitons. Section [5] summarizes the main conclusions of the paper. 



2. Theoretical model 



For the simple one- dimensional case, the optical field satisfies the Helmholtz 
equation 

^ ^ + kln\E)E = 0, (1) 



oz 2 ox 2 
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where k = oj/c is the propagation constant in vacuum, c is the speed of light, 
and u is the frequency. Introducing a normalization appropriate to a forward 
propagating beam, E(X,Z) = A(X, Z)exp(ikZ), a nonlocal nonlinearity 
n(E) = n + An and assuming that the approximation n 2 (E) = n 2 , + 2noAn, 
where k is the wave number in the medium, and An is the nonlinear induced 
change of the refractive index, which satisfies 

d 2 An 

™ 2 m^T-^ + n 2 \E\=0, (2) 

where n 2 is the nonlinear coefficient. We have employed the following nor- 
malizations, z = Z/L D , x = y/2X/wo, u(x,z) = ^ykn2Ld/n A(x } z), An = 
no4>/(kL£)), d = d/y/2. Wq is a transverse scale parameter that we shall 
later relate to the width of nonparaxial soliton beams. This scale parameter 
can also be considered as equivalent to the waist of a (reference) paraxial 
Gaussian beam, at z — 0, which has a diffraction length Ld = kw\j2. w m is 
the characteristic length of the nonlinear response, and d = w m /wQ stands for 
the degree of nonlocality of the nonlinear response. So, we get the following 
non-paraxial nonlinear Schrodinger equation(NNSE) for the dimensionless 
amplitude u of the light field coupled to the equation for normalized nonlin- 
ear induced change of the refractive index 



d 2 u du 1 d 2 u 
^ + Z ^ + 2^ + <M = °' (3a) 



^-0+M 2 = O, (3b) 

where k = 1/(/c 2 Wq) is the non-paraxial paremeter of the NNSE. In the limits 
k — y and d — > 0, the nonlinear Schrodinger equation(NSE) can be recovered 
from the system ([3]) which describes a local nonlinear response at d — > and 
a strongly nonlocal response at d — > oo. For Eq.( [3b|) . we can also write it 
into the form of convolution, 

/+oo ^ ^ 

R(x — x )I(x )dx , (4) 
-oo 

where I = I(x, z) = \u(x, z) \ 2 , R{x) = l/(V2d)exp(-y/2\x\/d). 
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3. Exact analytical soltions solutions 



When the nonlocal response is high, i.e., when the response function R(x) 
is wider compared with the extent of the beam. So, to find exact analytical 
solutions for the spatial solitons, we expand R(x , z) around the point x = x 
and combine Eqs. (l3a|) and (J4J to obtain j2l| 

d 2 u .du ld 2 u P 1 „ 2 

where \x = —R > [R' = d 2 R(x)\ x=0 , R < because R is a maximum of 
R(x)], P = f \u(x, z)\ 2 dx is the beam power, and P is the input power at 



0. We search for a solution to Eq. ([5]) of the Gaussian function form [21 



y r Poexp[ie{z)] x 2 . 2 

u{x, z) = ===== exp[- + lT W x ]' ( 6 ) 

/ 7TW{Z) ZW{Z) 



where 9 is the phase of the complex amplitude of the solution, w is beam 
width, t represents the phase-front curvature of the beam, and they are all al- 
lowed to vary with propagation distance z. The real amplitude of the solution 



has the form y/Po/ y \^ttw, owing to the conservation of the power. Insert- 
ing the trial function above into Eq. (j5]), we obtain that r = 0, w = Wq = 



l/Gu^Po 174 ), and 9 = (-dw 2 ± yf (Pvfi - 2d 2 w 2 K + 2y/2dwfaP ) / (2dw 2 K) . 
So, we get the exact solitons solutions. When d = 2.01, the solution is shown 
in Fig. [1] (a). To illustrate the stability of the soliton, we do the direct numer- 
ical simulations of Eq. flH]) with input condition u\ z= o = v(l + p) exp(— iVx), 
where p(x) is a broadband random perturbation, v is the nonlocal nonparax- 
ial soliton solution as an initial condition, and V denotes transverse velocity 
parameter. From Fig. [1] (b), we can easily see that the soliton is stable. 

When the nonlocal response is weak, i.e., when the response function 
R(x) is narrower compared with the extent of the beam. So, to find exact 
analytical solutions for the spatial solitons, we expand I(x , z) around the 
point x = x and combine Eqs. ( |3al) and (j4j) to obtain 



d 2 u .du ld 2 u .. l9 d 2 \u\ 2 . „ 

where 7 = 1/2 R{x)x 2 dx. For weakly nonlocal meda 7 1 is a small 
nonlocality parameter. We search for a stationary soliton solution to Eq. (J7j) 
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of the form u(x, z) = v(x) exp(i(3z), where the profile v(x) is real, symmetric, 



and exponentially localized and the propagation constant f3 is positive [22 
Substituting this solution into Eq. ([7j), we can get 

g + 2(v 2 + ^ - P )v + 2 1V ^ = 0, (8) 
Eq. (jSJ) can be finally integrated to give the equation [22 



1 (7 

±x = — tanh _1 (— ) + 2 v ^tan~ 1 (2 v /7a), (9) 

where Vq = 2((3 — k(3 2 ), a 2 = (vq — v 2 ) / (1 + i^v 2 ) . This implicit relation gives 
the profile of nonlocal nonparaxial spatial solitons propagating in weakly 
nonlocal Kerr-like media, when 7 = and k — 0, i.e., for the paraxial case 
with local nonlinearity, v(x) = v sech(vox) . The exact soliton solution is 
shown in Fig. [1] (c) when d = 0.11. Using Eq. (jSJ), we can analytically find 
the power 

P(fi) = v + L^M tow - 1 (2V7»o), (10) 

The derivative dP/dfi can be gained from Eq. (JTUl) easily and it transpires 
that the power monotonically increases with the propagation constant, in- 



dicating that the solitons are stable [22|, |23j. To further demonstrate the 



stability of the exact soliton, we also do the direct numerical simulations of 
Eq. (J7|) and find that the soliton is stable (Fig. Q] (d)). Here, v(x) is the exact 
soliton solution as an initial condition. 

4. Numerical solitons solutions 

To discuss the influence of the degree of nonlocality and nonparaxiality on 
the width and stability of solitons, we use split step Fourier method (SSFM) 



and spectral renormalization method [24( to obtain the nonlocal nonparaxial 



soliton solutions by solving the equations (j3j), where a solution guess for field 
distribution is u(x) = sech(x). Fig. [2] (a) shows the soliton solutions with 
the different degree of nonlocality d. We find that the solitons power will 
increase with the degree of nonlocality d. However, as shown in Fig. |2] (b), 
when d — 0.11, d — 1.01 and d = 2.01, the peak values of the solitons intensity 
Mmax versus propagation distance are almost straight horizontal lines, that is 
to say, d does not influence the stability of solitons. To illustrate the stability 
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of the solitons, we do the direct numerical simulations of Eqs. (J3]) with input 
condition u\ z= q = v(l + p) exp(— iVx), where v(x) is the numerically nonlocal 
nonparaxial soliton solution as an initial condition. The result in Fig. [2] (c) 
shows that the solitons can be stable when d = 0.11, which is consistent with 
the exact soliton solution. 

Next, we discuss the influence of non-paraxial parameter k on the solitons 
properties. From Fig. [3] (a), we can also see that non-paraxial parameter k 
has no effect on the field distribution when V and d are constant. However, 
k can influence the stability of solitons. In Fig. [3] (b), we easily find that 
the peak values of the solitons power |w|maz decrease with the propagation 
distance when k = 5x 10~ 3 and k = 1 x 10~ 2 . The solitons will be no more 
stable when k = 5 x 10 -3 and k = 1 x 10~ 2 . Also, we get the same conclusions 
by comparing Fig. [2] (c), Fig. [3] (c) and (d), that is, the solitons will be no 
more stable when the non-paraxial parameter k is bigger. When k is bigger, 
the width of nonparaxial soliton beams wq should be smaller. The bigger k 
can introduce the first type of nonparaxiality, so the solitons are no stable. 

5. Conclusion 

In conclusion, we investigate the nonparaxial solitons in the nonlocal 
nonlinear medium. We find an exact analytical solution to the nonlocal 
nonlinear nonparaxial propagation equation using Taylor expansion, when 
nonlocality is high and weak. Moreover, we numerically find that the degree 
of nonlocality and nonparaxiality have different influence on the width and 
stability of solitons. 
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List of Figures 



1 (Color online) (a) The profile of the exact soliton solution 
when d = 2.01. (b) Simulated propagation of the beam at 
d = 2.01 with p = 2% noise for the exact soliton solution 
initial condition, (c) The profile of the exact soliton solution 
when d = 0.11. (b) Simulated propagation of the beam at 
d = 0.11 with p = 2% noise for the exact soliton solution 
initial condition. The other parameters k — 1 x 10 -3 and 
V = 10 

2 (Color online) (a) The profiles of the solitons. (b) Evolution 
of the peak value of the solitons power H 2 ^. ( c ) Simulated 
propagation of the beam at d = 0.11 with p = 2% noise for a 
numerical soliton solution initial condition. The other param- 
eters k = 1 x 10~ 3 and V = 10 

3 (Color online) (a) The profiles of the solitons. (b) Evolution 
of the peak values of the solitons power |w|^ a:r . (c) and (d) 
show the simulated propagations of the beams with p = 2% 
noise at « = 5 x 10" 3 and k — 1 x 10~ 2 , respectively. The 
other parameters V — 10 and <i = 0.11 
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Figure 1: (Color online) (a) The profile of the exact soliton solution when d = 2.01. (b) 
Simulated propagation of the beam at d = 2.01 with p = 2% noise for the exact soliton 
solution initial condition, (c) The profile of the exact soliton solution when d = 0.11. (b) 
Simulated propagation of the beam at d — 0.11 with p = 2% noise for the exact soliton 
solution initial condition. The other parameters k = 1 x 10~ 3 and V = 10. 
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Figure 2: (Color online) (a) The profiles of the solitons. (b) Evolution of the peak value of 
the solitons power |u|m aa: . (c) Simulated propagation of the beam at d = 0.11 with p = 2% 
noise for a numerical soliton solution initial condition. The other parameters n = 1 x 10 -3 
and V = 10. 
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Figure 3: (Color online) (a) The profiles of the solitons. (b) Evolution of the peak values 
of the solitons power \u\ max- ( c ) an< ^ (d) show the simulated propagations of the beams 
with p — 2% noise at n = 5 x and k = 1 X 10~ 2 , respectively. The other parameters 
V = 10 and d = 0.11. 
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